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Let  H^(u}  s)  denote  the  harmonic  polynomial  of  degree  at 
most  n  found  by  Interpoletirai  in  2n  1  points  to  a  function 
u  given  on  the  boundary  C  of  a  region  D  of  the  conplex  emplane* 

It  is  proved  that  (a)  for  any  bounded  D  there  always  exist 
interpolation  points  on  C  so  that  can  be  uniquely  determined 
for  each  n,  and  (b)  for  a  wide  class  of  Jordan  regions  D  and 
for  boundary  data  u  with  a  smooth  first  derivative  on  C  the 
points  of  Interpolation  on  C  can  be  chosen  so  that  i^V*»  •) 
e^dstsy  z  e  Q-*'Dy  and  gives  the  solution  of  the 
Dlrlchlet  problem  for  u  and  D.  Explicit  formulas  are  derived 
for  in  the  case  of  interpolation  on  a  circle  and  on  an 
ellipse}  and  convergence  is  proved  in  these  cases  for  arbitrary 
continuous  boundary  data.  Various  generalizations  are  indicated. 
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1.  Introduction .  Let  U  be  a  function  on  an  open  set  D  in  the 
complex  s-plane,  t  -  x  iy,  to  the  real  numbers.  The  function  U 
Is  said  to  be  harmonic  on  D  If  It  is  continuous  there  together  with 
Its  partial  derivatives  of  the  first  two  orders  with  respect  to  x  and 
y,  and  If  It  satisfies  Laplace's  equation 


A  function  Is  said  to  be  harmonic  at  a  point  If  It  Is  harmonic  in 
some  neighborhood  of  the  point.  In  what  follows,  we  generally  shall 
taki;  D  to  be  a  region,  by  which  we  mean  a  non-ompty  open  connected 
subset  of  the  plane. 

Suppose  now  that  on  the  boundary  C  of  a  region  D  a  continuous 
function  u  from  C  to  the  real  numbers  Is  given.  It  is  well  known 
that  if  D  satisfies  certain  mild  restrictions  -  for  example.  If  its 
coiiq)lement  is  such  that  no  component  reduces  to  a  point  -  then  there 
exists  a  unique  function  U  harmonic  on  D  and  continuous  on  D  ♦  C 
which  coincides  with  u  on  C.  The  construction  of  this  function  U 
(or  In  purely  theoretical  contexts,  the  proof  of  Its  existence)  Is 
the  substance  of  the  famous  Dlrichlet  problem,  or  first  boundary 
value  problem  of  potential  theory,**  There  are  many  applications^ 

*  This  research  was  supported  by  the  United  States  Air  Force  through 
the  Air  Force  Office  of  Scientific  Research  of  the  Air  Research  and 
Development  Command,  under  Contract  No,  AF  1*9(638)  -  862, 

An  up-to-date  introduction  to  the  problem  will  be  found  in  [l.  Chapter  V], 


and  generalizations  of  the  two-dimensional  formulation  given  here* 

A  good  many  years  ago,  J,L»  Walsh  [13]  [lit]  raised  a  question  as 
to  whether  it  might  not  be  possible  to  represent  the  solution  of  the 
Dirichlet  problem  as  the  limit  of  a  sequence  of  hannonic  polynomiala 
found  by  Interpolation  to  the  boundary  data  u  on  C  in  suitably  chosen 
points.  Although  a  number  of  constructive  methods  of  finding  solu¬ 
tions  of  the  Dirichlet  problem  are  now  known,  the  suggestion  of  Walsh 
might  yield  a  particularly  simple  analytic  approximation  process,  if 
only  a  reasonably  general  convergence  theory  for  it  could  be  provided* 

It  would  appear  that  this  has  never  been  done  heretofore.  In  the 
present  paper,  it  is  established  that  there  exist  convergent  sequences 
of  harmonic  Interpolation  polynomials  for  a  wide  class  of  regions  and 
boundary  data,  but  the  last  word  on  the  general  subject  has  not  been 
said  by  any  means. 

In  the  references  given  above,  Walsh  indicates  how  to  obtain  an 
affirmative  answer  to  his  question  for  the  special  case  In  which  C  is 
a  circle*  Recently  [10]  he  has  supplemented  this  with  a  solution  for 
the  ellipse,  but  it  is  one  in  which  it  is  necessary  to  impose  a 
smoothness  condition  on  the  boundary  data  u  beyond  continuity*  Convej>- 
gence  is  then  achieved  on  the  ellipse  as  well  as  interior  to  it* 
Practical  considerations  and  experience  with  similar  but  easier  problems 
in  complex  polynomial  interpolation  suggest  that  the  primary  goal  of 
a  general  convergence  theory  in  the  present  problems  might  well  be  to 
establish  convergence  inside  the  region  under  minimal  hypotheses  on 
u  and  C,  without  regard  to  convergence  on  the  boundary,  where  the 
problem  is  more  delicate  and  where  in  the  applications  the  sought-for 
fiinction  is  known  anyhow. 

A  survey  of  the  general  problem  and  of  the  analogous  one  for 
complex  polynomial  interpolation  has  recently  been  published  by  the 
author  [U],  in  which  a  necessary  and  sufficient  condition  for  conver¬ 
gence  based  on  linear  operator  theory  is  given*  But  the  condition 

*  For  instance,  the  gap  in  theory  is  noted  by  Krylov  [17,  p.  U89].  He 
presents  in  [17]  a  formal  construction  of  the  Green's  function  by 
harmonic  polynomial  interpolation,  and  in  a  numerical  illustration 
it  seems  to  give  good  results. 
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gives  little  Indication  as  to  whether  In  a  given  Instance  there  exists 
a  convergent  sequence  of  harmonic  Interpolation  polynomials  and  as  to 
how  to  construct  It  If  so. 

In  the  present  paper  we  first  (Section  2)  examine  the  structure 
of  harmonic  interpolation  polynomials.  It  is  proved  that  for  any 
bounded  region  there  always  exists  a  point  set  S  on  the  boundary  such 
that  a  harmonic  polynomial  of  degree  at  most  n  which  assumes  2n  1 
preassigned  values  at  2n  ♦  1  of  these  points  S  will  be  uniquely 
determined  for  n>lj2,  ...  •  The  existence  question  requires  some 
attention  because  unlike  In  the  case  of  con^lex  polynomial  InterpolA-’ 
tlon  there  can  be  situations  involving  distinct  points  of  interpolation 
In  which  the  harmonic  Interpolation  polynomial  of  appropriate  degree 
Is  Indeterminate,  A  study  of  the  structure  of  the  Interpolation 
polynomial  Is  Important  because  the  basic  difficulty  In  the  whole  con¬ 
vergence  problem  for  harmonic  polynomial  Interpolation  Is  the  absence 
of  any  compact  formulas  like  the  Lagrange  and  the  Cauchy-Hermlte 
Interpolation  formulas  which  are  available  in  the  complex  polynomial 
case. 

In  Section  3  the  convergence  properties  of  harmonic  interpolation 
polynomials  are  related  to  other  types  of  harmonic  polynomial  approxln»- 
tlons.  It  is  thereby  proved  3»U)  that  for  any  region  with  a 

reasonably  smooth  boundary  C  and  for  any  boundary  data  u  possessing  a 
smooth  first-  derivative  it  is  possible  to  solve  the  Dirichlet  problem 
in  teims  of  the  limit  of  a  sequence  of  harmonic  polynomials  found  by 
interpolation  to  u  on  the  boundary.  The  convergence  takes  place  in 
the  closed  region,  Oenerallzatlons  and  situations  in  which  restrictions 
on  the  boundary  can  be  eliminated  are  indicated.  The  results  of  this 
section  do  give  an  affirmative  answer  to  Walsh's  original  question  for 
a  wide  class  of  regions  and  boundary  values ^  but  they  do  not  provide 
very  explicit  instructions  as  to  how  to  choose  successful  interpolation 
points  in  practice. 

Finally  in  Section  U  we  consider  a  special  choice  of  interpolation 
points  which  is  the  natural  analogue,  for  a  general  boundary  curve,  of 
points  equally  spaced  on  a  circle.  For  the  cases  of  the  circle  and  the 
ellipse,  explicit  formulas  are  derived  for  the  harmonic  interpolation 
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polynomials .  The  derivation  is  based  on  the  classical  Faber  polynomials* 
r.onvergence  is  proved  for  continuous  boundary  values.  Although  Walsh's 
results  [18]  anticipate  these  to  some  extent,  the  formula  for  the  poly¬ 
nomial  for  the  ellipse  is  new,  as  is  the  proof  of  convergence  inside 
the  ellipse  with  no  smoothness  conditions  on  the  boundary  data  beyond 
continuity. 

The  methods  are  largely  elementary,  in  the  sense  that  they  are 
based  on  material  usually  covered  in  good  first  courses  in  classical 
complex  variable  theory  and  linear  algebra.  It  is  only  in  the  proof 
of  Theorem  3.U  that  references  to  sophisticated  work  on  approximation 
theory  and  conformal  mappings  are  needed. 


This  is  a  paper  on  theory,  not  on  practice,  but  it  does  contain 
at  least  one  Indication  as  to  how  the  theory  might  be  put  to  work  in 
the  computation  laboratory.  The  all-important  issue  in  the  practical 
applications  will  certainly  be  the  correct  choice  of  the  interpolation 
points  on  the  boundary  C  of  the  region.  It  is  shown  below  in  Section 
3,^^eo£ems  3.2,  3.3,  and  3.U,  that  in  many,  if  not  all  cases,  the 
successful  choice  of  the  points  for  the  interpolation  polynomial  of 
n-th  degree  will  be  a  choice  which  maximizes  the  absolute  value  of  a 
certain  determinant  which  first  appears  in  the  display  (2.3)  below. 
(The  polynomials  Pj^*)  which  appear  in  (2»3)  are  any  conveniently 
chosen  complex  polynomials  in  z  with  degrees  coinciding  with  their 
subscripts.  •  The  choice  Pj[®)  ”  would  be  perml-ssible.)  If  it 
seemed  worthwhile  to  do  so,  then  for  a  curve  C  of  a  given  shape  it 
would  certainly  be  feasible  to  calculate  the  maximizing  points  for 
some  of  the  lower  values  of  n  and  put  them  in  a  library  ready  to  use 
for  any  boundary  data  u  which  might  be  presented.  If  this  is  all  that 
is  done,  then  the  interpolation  polynomial  of  degree  n  would  have  to 
be  determined  for  each  new  boundary  function  by  solving  a  certain  set 
of  linear  equations  in  2n  1  unknowns  for  the  coefficients  of  the 
polynomial.  But  in  Section  2  a  formula  is  derived  for  the  interpola¬ 
tion  polynomial  of  degree  n  which  expresses  it  as  a  linear  combination 
of  certain  harmonic  polynomials,  there  denoted  by  Bj^(b),  k  -  1,  ..., 
2n  +  1,  which  do  not  depend  on  the  boundary  data.  The  coefficients 
in  the  linear  combination  are  merely  the  values  of  the  boundary  data 


- 


at  the  chosen  2n 1  points  on  C«  It  would  be  possible  to  calculate 
the  harmonic  polynomials  once  and  for  all,  after  the  points  of 

Interpolation  on  C  have  been  chosen,  and  to  put  them  In  the  library 
along  with  the  Interpolation  points.  One  would  then  have  an  exceedingly 
simple  and  rapid  method  of  finding  an  analytic  approximation  to  the 
solution  of  any  Dlrlchlet  problem  that  might  be  presented  for  that 
particular  region. 


2.  The  structure  of  harmonic  Interpolation  polynomials.  A  harmonic 
polynomial  of  degree  n  is  an  expression  of  the  form 

n  .  . 

h(B)  -  o^  I  (o^.  r'^  cos  je  +  r  sin  j8), 

o  j.i  J  J 

2  -  r(cos  ©  •*■  1  sin  ©), 

where  o<^,  oc,,  ...,  o<  ,  ,  ...,  ?  are  real,  oC  and  ?  are 

0  1  n  1  n  n  n 

not  both  zero,  and  1  la  the  imaginary  unit.  This  can  be  written  as 

(2,1)  h(z)  -  a  ♦  1  (a.  z^  ♦  a.  z^), 

>1  •>  ^ 


where  a^  ■  o^^,  a^  -  (o<j  -  ifj)/2,  J  ■  1,  ...,  n,  and  the  bar 

denotes  conjugate  complex.  The  right  side  of  (2.1)  is  the  real  part 
of  a  polynomial  in  z  of  degree  n  with  coefficients  a^,  2a^,  ...,  2a^, 
It  is  also  the  Imaginary  part  of  a  polynomial  in  z  of  degree  n  with 
coefficients  ia^,  2iaj^,  .,,,  2ia^,  Thus  every  harmonic  polynomial  is 
the  real  part  of  a  polynomial  in  z  and  also  the  Imaginary  part  of  a 
polynomial  in  z.  Conversely,  the  real  and  imaginary  parts  of  any 
polynomial  in  z  are  harmonic  polynomials. 


The  real  part  of  a  polynomial  in  z  is  at  the  same  time  the  real  part 

of  a  related  polynomial  in  z.  The  coefficients  of  the  latter  polynomial 

are  the  complex  conjugates  of  the  coefficients  of  the  former  polynomial. 

From  this  it  follows  that  the  real  part  (and  also  the  imaginary  part)  of 

2  -  -2 

any  finite  linear  combination  of  the  monomials  1,  z,  z  ,  ...,  z,  z  ,  ... 
with  complex  coefficients  is  also  a  harmonic  polynomial. 
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2*1  .  If  a  finite  linear  combination  of  the  monomials  1*  z, 
«2  "  " 

c  ,  •••,  z,  z  ,  vanishes  for  all  *  on  the  boundary  of  a  bounded 
region  D,  then  the  coefficients  are  all  »ero« 


For  the  proof,  suppose  first  that  the  linear  combination  Is  a 
harmonic  polynomial  h(z),  and  let  It  be  written  out  In  the  form  (2.1) • 
By  the  maximum  principle  for  harmonic  functions  [1,  p.  179]  h(z)  must 
vanish  Identically  on  the  region  D.  Therefore  the  complex  polynomial  in 
z  of  which  h(z)  Is  the  real  part  must  reduce  to  zero  or  to  a  pure 
Imaginary  constant  [1,  p»  69]*  This  implies  in  turn  that  all  the  coef¬ 
ficients  In  (24I)  must  be  zeros. 

More  generally,  let  the  linear  combination  be 

L  -  Z  a-  z**  +  Z  b.  z**  , 

J-0  J  J-1  ^ 


where  of  course  some  of  the  coefficients  may  be  zero  and  a^  now  is  not 
necessarily  real.  Then* 

(2.2)  L  -  R  (  2  a.  z^  ♦  £  E  z'^  ]  . 

0  1 


Since  L  and  therefore  Rl  vanishes  on  the  boundary  of  D,  the  polynomial 

In  z  In  the  square  brackets  in  (2.2)  reduces  to  a  pure-imaginary  constant, 

and  a.  ♦  ^ ■  0,  j  •  1,  •••,  n.  Therefore 
J  J 

1  —  —1 

•  L  -  loc+  Z  (aj  z*'  -  aj  z*^  ),  oC  real. 


n  1 

This  expression  is  the  imaginary  part  of  the  polynomial  lo^-t^  Z  2aj  s'", 
multiplied  by  1,  so  It  is  a  harmonic  polynomial  multiplied  ^ 
by  1,  and  since  It  vanishes  on  the  boundary  of  D,  all  of  Its  coefficients 
must  be  zero.  This  cong)letes  the  proof  of  Theorem  2.1  . 


If  Pj^(z),  P2(*)»  ...»  ?„(*)  »re  complex  polynomials  In  z  of 


respective  degrees  1,  2,  ...,  n,  then  the  expression 
n 

b  ♦  Z  (b.  p.(z)  +  b.  p,(z)  )  with  b  real  is  also  a  harmonic 

o  J  j  j  j  0 

polynomial  of  degree  n.  With  a  rearrangement  of  terms  It  can  be  put 


Into  the  form  of  h(z)  above. 


The  symbol  means  "real  part  of” 
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If  such  a  polynomial  assumes  given  values  u^i  u^y  •••) 
respectively  at  points  s^i  •••«  *2ntl*  coefficients 


^0*  V  •••*  ^n*  ^1*  ^n 

algebraic  equations 
n 
b 


are  a  solution  of  the  2n  1  linear 


0  *  Pjf'h’  *  ^3  P?y  )  ■  %'  h  .  1.  2,  ...,  2n*l, 

We  shall  now  examine  the  existence  and  uniqueness  of  the  solutions  of 
such  sets  of  equations. 


Given  the  system  of  linear  algebraic  equations 
n 


>’(%>  ■  'o  *  jJ/'j  *  “j  pj^’  •  ”h  ' 

h  *  1^  2  f  9  •  •  f  2n  ^  1^ 


where  the  numbers  are  regarded  as  given  complex  numbers  and  the 
letters  Cj  and  dj  stand  for  unknowns,  a  sufficient  condition  for  the 
existence  of  a  unique  solution  is  the  non-singularity  of  the  matrix 


(2.3)  A  - 


1  p^(z^)  ...  p„(t^)  ••• 

: 

:  :  • _ 

’  Pl^®2n+l^  **’  Pn^*‘2n+1^  Pl^^2n+1^  •**  V^®2n7 


which  is  equivalent  to  the  noi>- vanishing  of  the  determinant  det  A,  Now 

A  can  be  transformed  by  elementary  column  transformations  into  a  certain 

matrix  A  which  is  the  specialization  of  A  in  which  p.(z)  is  replaced  by  z^, 
0  J 

A  necessary  and  sufficient  condition  for  the  singularity  of  A^,  and 
therefore  for  that  of  A,  is  that  there  exists  some  linear  combination  of 
the  monomials  1,  z,  z”,  z,  with  coefficients  not  all  zero, 
which  vanishes  in  all  of  the  points  z^^,  ...,  *2n+l* 

such  a  linear  combination  is  a  harmonic  polynomial  of  degree  not  greater 
than  n.  Conversely  any  harmonic  polynomial  of  degree  not  greater  than  n 
is  such  a  linear  combination.  The  upshot  of  this  is  that  the  necessary 
and  sufficient  condition  for  the  non-vanishing  of  det  A  is  that  the  points 
2^,  ...,  *2n+l  shall  not  all  lie  on  an  algebraic  curve  given  by  an 
equation  of  the  type  H(z)  -  0,  where  H  is  a  harmonic  polynomial  of  degree 
n  or  less.  We  henceforth  proceed  under  the  assumption  that  this  non- 
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singularity  condition  Is  fulfilled.  In  specific  Instances  Itn  validity 
must  be  checked.  It  Is  an  Implication  of  the  condition  that  the  points 
*2n+l  ”*^st  all  be  distinct,  for  if  they  were  not,  the  rank  of 
A  would  of  course  be  less  than  2n  '«■  1  and  the  colxums  would  be  linearly 
Independent*  We  call  the  condition  n-s  for  short  ("n-s"  for  "nonr- 
singularlty*' )  • 

Suppose  that  in  the  system  of  equations  ^  discussed 

above,  the  numbers  u^^  are  all  real.  Then  P(Zj^)  •  ,  h  ■  1,  ••*, 

2n  +  l,  which  means  that  the  coefficients  Cj  and  dj  satisfy  the  equations 
n 


2n  +1.  Condition  n-s  now  implies  that  c^  "  c^,  Cj  -  dj,  j  ■  1»  ...,  n,  so 

n  _ 

P(z)  is  necessarily  of  the  form  6(2)  ■  b  +  Z  (b.  p.(z)  +  b.  p.(z)  ), 

o  J  J  J  J 

b  real.  We  formalize  this  iri  a  theorem: 

0  . —  . -  ..■■■■  ■— LI.  ■l.M, 

Theorem  2.2  .  Given  complex  polynomials  p^(z),  P2(z)>  •••> 

z  of  respective  degrees  1,  2,  ...,  n;  real  numbers  u^^,  U2, 

u^^^j  and  complex  numbers  z^,  Z2>  ..»>  ^2n+l  >^8 

condition ;  there  exists  a  unique  linear  combination  of  the  polynomials 

n  _ 

p.(2)  of  the  form  c  +  Z  (c.  p.(z)  +  d.  p .(z) )  which  assumes  the  value 
'*jj  in  the  point  h«  1,  2n  +  l,  Furthermore  d^  -  Cj,  so  this 

is  a  harmonic  polynomial  B(z)  of  degree  at  most  n, 

.l^egrem  2,3  .  In_JJh^orgm  2*2  if  the  numbers  u^^,  u^,  ..*,  'i2n+l 
respective  values  assumed  in  the  points  z^,  82n+l’  Py  *  harmonic 

polynomial  h(z)  ■  c^  +  f  (c^  z*^  ♦  Cj  z*^)  of  degree  at  most  n,  then 
B(z)~h(z) . 


The  proof  consists  in  first  putting  B(z)  in  the  form 
n  1 

a  +  Z  (a.  z*^  +  a.  z'^)  and  then  using  the  non-singularity  of  A  to  show 
0  ^  J  J  o 

that  all  the  coefficients  of  B(z)  -  h(z)  must  vanish.  The  theorem 

implies  that  there  is  essentially  one  and  only  one  harmonic  interpolation 

polynomial  B(z)  under  the  liypotheses  of  Theorem  2.2  regardless  of  the 

choice  of  the  base  polynomials  Pj. 


Consider  now  the  matrix  A  In  (2.3)  and  modify  It  by  replacing  the 

k-th  row-vector  with  the  row-vector  p  ■  (1,  Pj^(*),  p]T^»  •••> 

p  (z}  ).  Let  A. (z)  denote  the  new  matrix  so  obtained.  The  quotient 
n  K 

Bj^(z)  ■  det  A|^(s)  /  det  A  is  clearly  a  linear  combination  of  the  elentota 
of  ^  which  vanishes  for  z  -  Zj^,  Zg,  *k+l*  ****  ®2n+l  “*** 

equals  one  for  z  •  Zj^.  Therefore  by  Theorem  2.2  and  2.3  it  is  the 
unique  harmonic  polynomial  of  degree  at  most  n  which  behaves  in  this  way. 

It  follows  that  given  any  real  numbers  u^^,  u^,  *..>  U2n+1 
n  complex  polynomials  in  z  in  which  one  of  each  degree  from  1  to  n  is 


present,  the  unique  harmonic  oolynomlal  of  degree  at  most  n  which  assuaeo 


these  values  in  points  z^^,  ,,,,  *2n+l  condition  can 

be  represented  by  the  formula 


(2.U) 


H  (z)  -  Z  u.  B.  (z) 

n  ic-l  I' 


2n+l 
■  Z  u. 
k-1 


det  Aj^(z) 


If  we  write 

®k<*>  •  ‘’ko  *  J>'’k3  Pj**’  *  ^  > 

then  bj^^  is  the  cofactor  of  the  first  element  in  the  k— th  row  of  Aj^(z) 

(or  A)  divided  by  det  A,  bj^j,  j  -  1,  2,  n,  is  the  cofactor  of 

Pj(z)  in  Aj^(z)  divided  by  det  A,  and  J  **  is  tho 

cofactor  of  p .(z)  divided  by  det  A.  These  cofactors  divided  by  det  A 
J  —1 

form  the  k— th  column  of  the  matrix  A  ,  so  if  we  denote  this  coluim- 
vector  by  [A“^]j^,  we  can  write  Bj^(z)  -  p*[A”^]j^,  Thus  we  get  the 
compact  formula 

-1  -1 

H  (a)  -  Z  u.  p*[A  p  A  ^  (i, 

n  k-1 

where  Q  is  the  column-vector  with  elements  u^,  Ug*  •••>  '^n+1* 

We  conclude  with  a  result  which  is  basic  to  certain  convergence 
theorems  in  the  next  section. 


Theorem  2,U  .  Given  any  bounded  region  D,  for  each  n,  n  -  1,  2,  ,,,, 
there  always  exists  at  least  one  set  of  2n  + 1  points  on  the  boundary  of 


D  which  satisfies  the  n-s  condition. 
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The  statement  is  equivalent  to  saying  that  given  any  bounded  region 
D,  points  •••(  ®2n+l  exist  on  the  boundary  such  that  the 

matrix  in  which  the  k-th  row  is  (1,  z^)  f 

k  ■  1,  2n+l,  is  non-singular.  We  use  what  is  known  in  linear 

computation  theory  as  an  "escalator  method"  to  establish  this. 


Consider  first  the  upper  left-hand  (n  *  1)— rowed  principal  minor 
of  A^,  which  we  denote  by  V.  This  is  merely  a  Vandermonde  matrix  and  Is 
non-singular  whenever  the  points  are  distinct,  as  may  be 

shown  in  various  ways  -  for  example,  by  direct  evaluation  of  the  determinant. 
We  choose  any  set  of  distinct  points  on  the  boundary  of  D  as  our  points 
z^,  and  pass  to  the  consideration  of  the  minor 


V 


The  determinant  of  this  minor,  considered  as  a  function  of  ^ 

linear  confcination  of  the  monomials  1,  »^^2»  ®n+2’  *n+2*  ^  which 

the  coefficients  are  the  cofactors  of  the  elements  in  the  last  row.  If 
for  all  z^^2  ®n  the  boundary  of  D  this  determinant  were  to  vanish 
Identically,  then  bv  Theorem  2,1  the  coefficients  would  all  have  to  he 
zero.  But  the  coefficient  of  Zj^^2  zeroj  it  is  det  V,  Therefore 

there  must  exist  locations  for  z^^g  o'™  boundary  of  D  such  that 
det  ^Q^(®p+2^  ^  choose  one  such,  and  pass  to  the  next  principal 

minor  of  A^,  say  The  determinant  of  this  is  a  linear  combination 

V3*  V3»  ®n+3»  V3*  ooefficient  of  z^^^  is 

det  since  this  is  not  zero,  there  must  be  a  way  of 

choosing  o'™  the  boundary  of  D  so  that  det  ^2^®n+3^  /  0,  And  so 

on}  induction  completes  the  proof. 

^eorem  2  can  also  be  derived  from  a  more  general  result  outlined 
in  [7,  pp.  U87-I488].  This  is  to  the  effect  that  if  u^,  u^,  ...,  u^  are 
linearly  independent  real— valued  functions,  and  if  for  any  n  points  on  an 
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arc  oC  it  is  always  possible  to  pass  a  "curve”  with  equation  ♦ 

...  -t'  a  u  (a)  ■  0  through  the  points,  then  o<  must  coincide  with  such 
a  curve .  In  our  case ,  if  no  sot  of  2n  + 1  points  satisfying  the  n~B 
condition  were  to  exist  on  the  boundary  of  D,  then  each  such  point  set 
must  lie  on  a  locus  H(e)  ■  0,  where  H  is  a  harmonic  polynomial  of 
degree  at  most  n  with  coefficients  not  all  zero.  But  then  according  to 
the  above  result  the  boundary  of  D  would  have  to  coincide  with  such  a 
curve,  and  this  would  be  impossible  by  the  maximum  principle  for  harmonic 
functions . 

3.  Convergence  theorems  for  sequences  of  harmonic  interpolation 
polynomials «  Let  D  be  an  arbitrary  botinded  region  and  let  C  be  its 
boundary.  Consider  the  following  infinite  sequence  of  point  sets  lying 
on  Ct 

^1  ■  "ij’  ‘13’ 

Sj  ■  **n’  *22'  *23’  *21l’  *25* 

®n  “  ^*nl»  ®n2»  *n,2n4l^ 

We  assume  throughout  this  section  that  for  each  n,  satisfies  the  n-e 
condition  of  Section  2.  Let  k  ■  1,  ...,  2n>l,  denote  for 

each  k  the  (iinique)  harmonic  polynomial  of  degree  at  most  n  which  vanishes 
at  all  points 'of  except  the  point  at  which  it  equals  one.  For 

any  function  u  given  on  C,  we  construct  the  harmonic  polynomial  of  degree 
at  most  n  found  by  interpolation  to  u  at  the  points  S^j  written  in  the 
form  (2.U)  this  polynomial  is 

2n+l 

(3.1)  i)  .  Z  u(.  )  B  (z), 

k"l 

The  method  used  in  this  section  to  study  the  convergence  properties 
of  the  sequence  H^,  H^,  ...,  consists  in  referring  them  to  the 
properties  of  other  sorts  of  harmonic  polynomial  approximations  to  u, 

V/e  recall  that  if  h  is  any  harmonic  polynomial  of  degree  not  more  than 
n,  then  by  Theorem  2.3, 

2n+l 

H^(h}z)  .  2  h(V^ 

k"l 
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Suppose  now  that  at  least  for  a  particular  u  it  is  possible  to  solve  the 
Dirichlet  Problem  for  u  and  D|  let  U  denote  the  solution.  For  any  z 
in  the  finite  plane  for  which  U  exists, 

2n+l 

K^(ui  z)  -  U(*)  - 

2n+l 

Now  let  T  (z)  ■  I  |B  ,,(*)!•  By  the  maximum  principle  for  harmonic 
n  nx 

functions , 

|h(z)-U(z)I  <  |u(z)  -  h(z)|,  zeD+C. 

This  with  (3.1)  implies  that 

(3.2)  |H„(u;  *)  -  U(*)l  <  [  z  e  D  ♦  C. 

We  base  our  convergence  theorems  all  on  this  inequality.  In  the 
language  of  linear  operator  theory,  the  function  T  (z)  is  the  total 
variation  of  the  kernel  of  H^,  regarded  as  a  linear  operator  j  or  again, 
it  is  the  norm  of  the  functional  H  on  the  space  of  functions  u»  Clearly 
it  must  play  a  key  role  in  the  convergence  theory. 

In  what  follows  we  shall  have  frequent  occasion  to  make  use  of  the 
"o”  and  "0"  notation  when  comparing  orders  of  magnitudes  of  functions 
defined  on  the  positive  real  integers.  The  reader  will  doubtless  recall 
that  the  notation  a(z,  n)  -  o(b(n))  means  that  ^Ijj^  a(z,  n)/^(n)  ■  0, 
Also,  a(z,  n)’  ■  0(b(n))  means  that  a  constant  M  exists  which  is  independent 
of  n  such  that  a(z,  n)  <  M  b(n)  for  all  n  sufficiently  large.  A 
uniformity  condition  with  respect  to  z  on  a  set  Z  means  that  the  limit 
in  the  '*0”  case  is  uniform  with  respect  to  z  and  the  M  in  the  "0"  case 
is  independent  of  z  for  z  on  Z, 

Theorem  3.1.  Let  C  be  a  Jordan  curve  ,  let  u  be  continuous  on  C,  and  let 
T^(z)  ■  0(1)  uniformly  on  the  subset  Y  of  D+C.  Then  ^Ijj^  "  U(z) 

«  For  a  further  development  of  this  point  of  view,  see  [hj.  The  comparison 
technique  in  (3.2)  is  a  familiar  one  in  the  theory  of  interpolation j  for 
example  see  [12,  Chapter  XIV],  [6,  Chapter  IV], 

A  Jordan  curve  (or  simple  closed  curve)  is  homeomorphic  to  a  circle. 
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uniformly  on  V  .  If  is  unbounded  at  any  point  then  there 

exists  a  continuous  function  u  for  which  {H  (u:  z  )}  is  unbounded* 

■  ■  — '  ■  '  "  n  o  . —  — 

The  theorem  was  first  announced  in  [U].  The  convergence  statement 
follows  from  (3*2)  and  a  theorem  of  Walsh  [17,  p*  169]  to  the  effect  that 
given  any  continuous  function  u  on  a  Jordan  curve,  there  exists  a  sequence 
of  harmonic  polynomials  h^(z),  h2(z),  ...  of  respective  degrees  not 
greater  than  n  which  converges  to  u  uniformly  on  C,  The  last  sentence 
follows  from  linear  operator  theory;  see  [U]. 

Corollary.  K  "  0(l/n),  k  -  1,  ...,  2n+  1,  uniformly  in 

k  and  uniformly  for  z  on  any  closed  subset  of  the  Jordan  region  D,  and 
if  u  is  continuous  on  C,  then  11^^  exists  imiformly  for  z 

on  any  closed  subset  of  D  and  provides  there  the  solution  of  the 
Dirichlet  Problem  for  u  and  D, 

The  hypothesis  of  course  implies  that  T^(z)  is  uniformly  bounded 
on  any  closed  subset  of  D.  We  have  stated  the  corollary  formally 
because  a  study  of  various  special  cases,  such  as  the  ones  discussed  in 
Section  li  below,  seems  to  indicate  that  on  any  sufficiently  smooth 
Jordan  curve  there  always  exists  a  sequence  of  point  sets  S^^,  Sg,  ... 
such  that  the  condition  on  in  the  hypothesis  of  the  corollary  is 
satisfied • 

It  is  to  be  noted  that  under  the  hypotheses  of  Theorem  3*1,  the 
degree  of  convergence  to  u  on  D  of  any  comjjarison  sequence  h^(z), 
h2(z),  ...  is  duly  reproduced  by  Hj^(u;  z),  H2.(u}  z),  ...  on  y  . 

For  example,  if  a  sequence  h^(z),  h2(z),  ...  exists  such  that 

u(z)  -  h^(z)  ■  0(n“^) ,  q  >  0,  uniformly  on  C,  then  H^(u;  z)  —  U(z)  -  0(n“^) 

uniformly  on  y  ,  It  is  known  [3],  [19],  that  if  C  is  analytic*  and  if  u 

(k) 

has  a  k— th  derivative  u  with  respect  to  arc  length  on  C  which  satisfies 
the  following  condition  (called  a  Lipschitz  condition  of  order  oC) » 


*  For  any  Jordan  curve  C  there  exist  parametric  representations  of  the 
form  z  -  f,(9)  +  i  f2(®)»  where  fj^  and  f2  are  real  continuous  functions 
of  period  2tT,  and  where  any  two  solutions  0  of  the  equation 
for  a  given  z  on  C  differ  by  an  integral  multiple  of  0.  A  Jordan  curve 
is  said  to  be  analytic  is  there  exists  such  a  parametric  representation 
in  which  f^  and  f2  are  analytic  functions  of  0  and  |f^(0)|  ♦  |f2(®)|  /  0* 
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X  >  0, 

0  <  oC  <  1,  all  and  Zg  on  C, 
then  there  exist  harmonic  polynomials  exhibiting  the  degree  of  convergence 

The  problem  of  finding  point  sets  on  a  general  Jordan  curve  such 

that  T  (z)  is  uniformly  bounded  for  z  inside  the  curve  remains  open  as 
n 

this  is  being  written,  but  it  is  easy  to  show  that  there  exist  sets 
for  which  the  sequence  (T^(z)}  increases  no  more  rapidly  than  does 
2n  +  1  for  z  on  D  +  C ,  We  now  do  this • 

Theorem  3.2  .  l£t  P  be  any  closed  bounded  point  set  with  the  property 
that  for  each  n  a  subset  of  2n  ♦  1  points  satisfying  the  n-s  condition 
exists.  Then  there  exists  a  sequence  of  subsets  of  P  ,  S2,  ...  t 

each  satisfying  the  n-s  condition,  such  that  when  is  constructed 

for  the  n-th  subset,  n"l,  2,  ...  ,  zeP. 

For  the  proof,  consider  the  matrix  A  appearing  in  (2.3)  with 
Zj^,  ...,  *2n+l  ^  polynomials  p^,  ...»  in  A  may  be 

chosen  arbitrarily  in  what  follows,  except  that  as  usual  their  degrees 
must  agree  with  their  subscripts.)  Now  |det  A|  is  a  continuous  function 
of  the  2n  +  1  independent  variables  z^^,  ...,  *2n'*-l*  domain 

is  a  compaat  set.  Therefore  the  function  has  a  maximum  on  this  set. 

The  maximum  cannot  be  zero,  because  the  n-s  condition  in  the  Lemma  is 
equivalent  to  stating  that  det  A  /  0  for  at  least  one  S^.  Let 
be  any  one  of  the  sets  S  which  maximize  |det  A |  and  denote  the  corres- 
ponding  matrix  by  A  .  Consider  now  the  matrix  Aj^(z)  obtained  from  A 
by  replacing  the  k-th  row-vector  with  (1,  Pj^(z),  ...,  Pj,(*)»  P^T*7»  ...» 

) ) ,  Because  of  the  maximizing  property  of  S^,  it  must  be  true  that 
|det  A*(z)  I  <  |det  A*]  for  all  z  on  P  ,  Let  B^j^(z)  in  the  definition 
of  T^(®)  specialized  to  (det  A*(z))/(det  A*).  Then  i  ^ 

for  all  z  on  P  ,  and  the  inequality  in  the  Theorem  follows  from  this. 
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The  reasoning  above  is  reminiscent  of  studies  of  interpolation 
points  with  extremal  properties  which  have  been  made  by  Fekete,  Leja, 
Shen,  and  others  in  connection  with  complex  polynomial  interpolation; 
see  (17,  66  7«V— 7«8]« 

Obviously  the  points  satisfy  the  n-s  condition  and  are  distinct. 

From  Theorem  2.U  it  can  be  seen  that  at  least  whenever  V  contains 
the  bom  iry  of  a  bounded  region,  the  hypotheses  of  Theorem  3.2  on  H 
will  be  fulfilled. 

Theorem  3.3  •  Let  D  be  a  bounded  simply  connected  region  for  which  the 
Dlrichlet  Problem  can  be  solved,  at  least  for  a  certain  f\inction  u 
continuous  of  the  boundary  C  of  D»  If  there  exist  harmonic  polynomials 
h^(z),  h2(z),  ...  of  respective  degrees  no  greater  than  1,  2,  ...  , 
such  that  u(z)  -  h^(z)  -  o(l/n)  uniformly  on  C,  then  there  exists  a 
sequence  of  point  sets  S^,  S 
exists  uniformly  on  D+  C  and  provides  the  solution  of  the  Dirichlet 
Problem  for  u  and  D. 

Theorem  2.U  guarantees  that  the  boundary  C  in  this  theorem  satisfies 

the  hypotheses  imposed  on  P  in  Theorem  3.2,  so  there  exists  a  certain 

sequence  of  subsets  of  C,  say  S*,  S^,  ...,  such  that  <  2n+l, 

z  e  D+C  when  T  (z)  is  constructed  for  S*.  Substituting  this 
n  n 

inequality  into  (3.2),  we  obtain 

(3.3)  (u;  z)  -  U(z) I  <  |u(z)  -  h  (z) |  (2n  ♦  2),  z  s  D  ♦  C, 

from  which  the  conclusion  of  Theorem  3.3  follows  at  once. 

The  polynomials  H  (uj  z)  here  interpolate  to  the  boundary  data  u 

^  I  . 

in  the  special  points  which  maximize  [det  A|.  For  a  brief 

discussion  of  computational  aspects  we  refer  the  reader  to  the  last 
paragraph  of  the  Int.J-'ductlon,  Section  1  above. 

Under  what  conditions  on  D+  C  and  on  the  boundary  data  u  do 
harmonic  polynomials  exliihlting  the  degree  of  convergence  o(l/n) 

exist?  We  have  indicated  above  that  if  C  is  an  analytic  Jordan  curve 
and  if  the  first  derivative  of  u  on  C  satisfies  a  Lipschitz  condition  of 
order  o<,  then  the  desired  harmonic  approximation  is  avaiLable.  It  is 
possible  however  to  relax  the  requirement  of  analyticity.  By  reference 


2*  ...  on  C  such  that  l^^H^(u;  z 
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to  some  work  of  Sewell  [10][ll]  we  shall  now  prove  that  the  Dirichlet 
Problem  can  be  solved  by  harmonic  polynomial  interpolation  for  a  rather 
wide  class  of  Jordan  regions  when  the  boundary  data  belong  to  a  Lipschitz 

class* 

Let  C  be  a  Jordan  curve  and  let  K  be  the  region  of  the  extended 
plane  exterior  to  C.  We  consider  a  conformal  mapping  of  K  onto  a  region 
in  another  complex  plane,  which  we  call  the  w-plane.  There  exists  a 
function  z  -  <!•('')»  analytic  and  univalent  on  |w|  >  1,  which  maps 
|w|  >  1  conformally  onto  K  so  that  the  point  at  infinity  in  the  w-plane 
goes  into  the  point  at  infinity  in  the  z— plane •  By  the  Osgoodr-Caratheodory 
Theorem  [2,  p.  86]  this  function  can  be  extended  continuously  onto 
|w|  -  1  so  that  it  gives  a  topological  mapping  of  |w|  ■  1  onto  C, 

Definition.  The  curve  C  is  of  Type  W  provided  that  l))  (w)  /  0  on 
Iw]  -  1  and  is  continuous  on  |w|  ■  1, 


The  primes  indicate  differentiation.  Here  and  in  what  follows,  the 
derivative  of  a  function  f  given  on  a  Jordan  curve  C  means  the  derivative 
in  a  one-dimensional  sense  with  respect  to  z  on  C» 


df 

Hi 


fU^) 


z»z 


1 


lim 


f(z)  -  f(z^) 
z  "z , 


^(k).  ^ 


^1) 

dz 


k"2,  3,  •••>  ^^^1 


Theorem  3*U  •  Let  D  be  a  region  of  which  the  boundary  is  a  Jordan  curve 
C  of  Type  W.  Let  u  and  its  Tirst  derivative  u*  be  continuous  on  C  and 
let  u'  satisfy  a  Lipschitz  condition; 


lu'(zj^)  -  uVz^)  I  <  "X  |Zj^  -  Zg  1°^, 

A  >  0,  all  z^  and  z^  on  C, 

of  some  order  o< ,  0  <  oc  <  1.  Then  there  exists  a  sequence  of  point  sets 

Si,  Sg,  •«.,  onC  such  that  j^ip  ^  exists  uniformly  on  D  +  C 

and  provides  the  solution  of  the  Dirichlet  Problem  for  u  and  D, 

It  is  the  plan  of  the  proof  to  show  that  under  the  hypotheses  of  the 
theorem  there  exist  harmonic  polynomj.als  h^Cz),  h^(z),  ,,,,  of  respective 
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degrees  at  most  2,  3>  •••  such  that 


(3.li)  u(*)-h^(z)  -  0[(  ], 

uniformly  for  z  on  C.  Now 

/  log  n  xoC+l  _  1  (log  _  /lx 

^  ir“  ^  n 

n 

so  if  we  can  obtain  (3«h),  then  the  hypotheses  of  Theorem  3»3  will  be 
satisfied.  Theorem  3,h  will  then  be  an  immediate  consequence  of 
Theorem  3 •3, 


The  theorem  of  Sewell  [11]  whieh  we  need  is  the  following: 

Let  the  point  set  P  consist  of  a  finite  number  of  bounded  closed 
regions  bounded  by  the  set  of  mutually  exterior  Jordan  curves 
,.*i  each  of  'l.vpe  W,  For  each  j,  j  «  1,  m,  let  the  function 
fj  be  analytic  interior  to  Cj,  be  continuous  in  the  corresponding 
closed  region ^  and  possess  on  Cj  a  continuous  k-th  derivative  which 
satisfies  a  Lipschitz  condition  with  exponent  o<,  0  <  o<  <  1.  Then 

for  each  integer  n,  n  >  2,  there  exists  a  complex  polynomial 
of  degree  at  most  n  ^  z,  such  that 


(3.5)  f(z)  -  P^(z)  -  0[  (  ], 

uniformly  for  z  on  P , 


In  the.  present  application,  we  use  the  case  in  which  m  ■  1, 
k  -  1,  oC  <  1. 

We  need  the  following  facts  concerning  C,  which  are  implied  by  the 
definition  of  Type  W : 

The  curTC  C  has  a  tangent  at  each  point.  If  T"(s)  denotes  the 
tangent  angle  with  the  positive  real  axis,  expressed  as  a  function  of 
arc  length  s  on  C,  then  T^(s)  satisfies  a  Lipschitz  condition  with 
exponent  unity  uniformly  in  s,  (See  [11]  for  further  explanation,)  Let 
z  -  (J(w)  be  analytic  and  univalent  for  |w|  <1  and  map  |w  |  <  1 
conformally  onto  D;  further  let  z  -  (^(w)  be  extended  continuously 
onto  |w(  ■  1.  Then  z  ■  (j)(e  ),  0  <  S  <  2it,  is  a  parametric  equation 
for  C«  Finally  let  w  ■  X(z)  be  the  inverse  of  z  ■  (J(w),  By  a 
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theorem  of  Kellogg  [8][5,  pp»  3h-35]  the  Lipschitz  condition  on  T^s) 
implies  that  ()i*(w)  and  "X-'C®)  exist,  and  that  for  each  oC,  0  <  <  1, 

(jl'(w)  satisfies  a  Lipschitz  condition  in  w  on  |w|  -  1  with  exponent  '><, 
and  xV*)  satisfies  a  Lipschitz  condition  in  z  on  C  with  exponent  cxT. 

In  particular,  <))'  is  continuous  on  |w|  ■  1  and  Tc*  is  continuous  on 
C,  From  the  continuity  of  the  derivatives  it  follows  that  positive 
constants  a  and  b  exist  such  that 


(3.6) 


(3.7) 


-  (^(Wg) 


<  a, 


all  |w^|  -  1,  Iwgl  -  1, 


<  b,  all  z^,  Zg  e  C, 


(The  existence  of  such  bounds  is  easily  proved  by  an  indirect  argument. 
For  a  detailed  direct  proof  which  goes  into  the  structure  of  the  bounds, 
see  (19,  pp.  38^386].) 


We  return  now  to  the  function  u  of  Theorem  3.U.  The  mapping 
z  ■  ()(w)  carries  u(z),  z  on  C,  into  u((J(w))  ■  v(w),  |w|  ■  1.  Now 


v'(w) 


du  du  dz 
clw  "  ^  *  3w 


u’(z)  (>'(w),  |w|  -  1, 


and 

v'(w^)  -  v'^Wg)  -  (u'(Zj^)  -  u*(z2))  (^'(w^)  ♦  u'(z2)((^'(Wj^)  -  (^'(Wg)), 

®1  "  *2  ■■  ‘^(''2^’  *1’  *2  ® 


liSt  be  an, upper  bound  for  |()l'(w)l,  |w|  ■  1,  and  let  Mg  be  an  upper 

bound  for  the  modulus  of  the  continuous  function  u*  on  C.  Also  let 
be  the  constant  in  the  Lipschitz  condition  satisfied  by  (J'(w)  on  |w|  ■  1, 
Then  we  obtain: 

(3.8)  Iv'(Wj^)  -  ▼'(Wg)  I  <  |z^-Zg|°^+  Xj^Mg  Iwj^-Wgl*' 


-  Xm. 


()»(w^)  -  (|>(Wg) 


Wi  -  Wg 


oC 

K  ~  K"  '^2 1°^ 


<  (A.M^  a°*^  +  “  ^^2  0  <  0^  <  1 
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(For  the  last  Inequality  we  used  (3»6).)  Thus  v'(w)  satisfies  a 
Lipschlts  condition  of  order  o<  on  |v| 

A  theorem  of  Walsh,  Sewell,  and  Elliott  [19,  p.  388],  based  on  work 
of  Privalov  [20,  vol,  I,  pp»  121-122]  on  conjugate  functions,  now 
guarantees  the  existence  of  a  function  g  analytic  in  the  complex  variable 
w  for  Iwl  <1,  continuous  on  [wj  <  1,  and  with  the  properties  that 
its  real  part  coincides  on  lw(  ■  1  with  v  and  also  that  its  derivative 
g'  satisfies  a  Lipschitz  condition  of  order  oC  on  jw)  -  1»  The  mapping 
w  -OcCz)  carries  g(w)  into  g('5C(z))  ■  f(z)  which  is  analytic  for  s 
on  D  and  continuous  on  D  +  C.  Also,  5^f(z)  ■  u(z),  z  on  C,  The 
derivative  f'  exists  at  each  point  of  C,  because 

-  g'(w^)  ■>c'(zj), 

w-  X(z),  Wj^  -  s,  e  C, 

Finally,  as  in  (3.8), 

If'(z^)  -  f'(z2)|  -  |(g'(w^)  -  g’(w2))  X'(z^)  ^  g'(w2)(X.'(z^)  -X'(Z2))| 

<  ^2^3  ^'^l  ~  ^^2  ^  ^  ^*1  ~  *2  ^ 

(X2M3  b*^  ^  *2!°^* 

-  X(Z;,^),  W2  -  X(Z2),  Iwj^l  -  IW2I  -  1, 

where  and  are  the  Ilpschitz  constants  of  g'  and  "X.'  respectively, 
and  are  upper  bounds  for  [X*  I  and  |g'|  respectively,  and  b  is  the 
bound  which  appears  in  (3.7). 

This  function  f  therefore  satisfies  the  condition  of  Sewell's  theorem 
with  m  ■  1,  k  ■  1,  c<  <  1.  Therefore  there  exists  a  sequence  of  conplex 
polynomials  in  z  which  satisfy  (3.5).  Let  h^(z)  in  (3.lj)  be  HPjj(z) 
in  (3.5).  We  have,  since  u  ■  Rf, 

|u(z)-h^(z)|  <  Ir(z)  -  P^(z)  I, 


f(z)  -  f(Zj^) 
z  -  z. 


g(w)  -  g(w^)  X(»)  -  X(Zj^) 


w  —  w, 


z  -  z. 
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and  the  existence  of  the  required  harmonic  polynomial  approximation  to  u 
follows  from  this  inequality.  This  concludes  the  proof  of  Theorem 

The  generality  of  the  point  set  P  in  Sewell's  result  and  in  our 
Theorem  3*2  suggests  that  Theorems  3*3  and  3«U  remain  true  if  D  is  replaced 
by  a  point  set  consisting  of  a  finite  number  of  mutually  exterior  bounded 
regions,  D^,  D2,  D^,  with  functions  u^^,  U2,  u^^  given 

respectively  on  their  boundaries.  (In  the  case  of  Theorem  3»U,  the  region 
Dj^  would  be  bounded  by  curves  of  Type  W.)  This  is  in  fact  true.  The 
limit  of  the  sequence  of  harmonic  interpoLation  polynomials  then  simul— 
taineously  provides  the  solution  of  the  Dirichlet  Problems  for  and 
with  a  on  Dg  and  U2  with  z  on  Dg,  and  so  on. 

The  result  of  Sewell  used  to  establish  (3.U)  and  thereby  Theorem  3»U 
seems  to  be  nearly  the  best  possible,  in  the  sense  that  under  the  stated 
conditions  on  the  approximated  function  f,  the  conditions  on  the  curve  C 
are  nearly  minimal  for  achieving  the  exhibited  degree  of  approximation. 
However,  the  hypotheses  of  Theorems  3.3  and  3oU  are  strong  enough  to 
guarantee  convergence  on  the  boundary  C  as  well  as  on  D,  and  it  was  pointed 
out  in  the  Introduction  that  because  of  this,  in  a  sense  these  theorems 
overshoot  the  goal  of  the  general  convergence  theory.  We  here  amplify 
the  remarks  in  the  Introduction  concerning  the  extra  difficulties  on  the 
boundary. 

Consider  the  case  in  which  C  is  the  unit  circle,  A  glance  at  the 
first  formula  in  Section  2  above  will  reveal  that  in  that  case,  harmonic 
polynomial  interpolation  on  C  reduces  to  interpolation  with  trigonometric 
sums.  In  this  case  it  is  known  [20,  vol,  II,  p.  37],  [6,  p,  120]  that  a 
sequence  of  sets  S*,  S*,  ...  exists  such  that  "  0(log  n) ,  [z]  ■  1, 

which  is  much  more  favorable  to  convergence  than  the  0(n)  given  to  us  by 
Theorem  3.2.  The  points  in  are  equally  spaced  on  |z|  ••  1.  But  to 
obtain  convergence,  additional  smoothness  conditions  are  needed  for  the 

*  Mergelyan  [9,  p.  8ljJ,  using  hypotheses  on  C  which  differ  slightly  from 
those  of  Sewell,  arrives  at  a  slightly  weaker  degree  of  convergence,  and 
shows  that  his  result  is  the  best  possible  for  the  class  of  functions  f 
and  class  of  curves  which  he  considers. 
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boundary  values  beyond  mere  continuity,  for  it  is  also  known  [20,  vol»  II, 

pp.  hh  ff .]  that  if  the  point  z  »  1  belongs  to  each  S^,  then  there 

exists  a  function  u  continuous  on  |z|  «»  1  such  that  {H  (uj  z)}  diverges 

for  all  z  /  1,  where  H^(u;  z)  is  found  by  interpolation  to  u  in  S^, 

For  such  a  function  u  the  Fourier  series  may  converge  uniformly  (20,  vol.  II, 

p.  b7],  allhotigh  (speaking  qualitatively)  the  partial  sums  of  the  P'ourier 

series  generally  do  not  give  a  particularly  good  trigonometric  polynomial 

approximation  to  a  merely  continuous  function.  We  shall  show  in  the  next 

section  that  for  any  continuous  u,  the  above  sequence  [H  (uj  z)}  does 

n 

converge  on  |z|  <  1,  no  matter  how  badly  it  behaves  on  lz|  “  lo 

In  defense  of  Theorems  3*3  and  3»U  it  might  be  said  that  there  are 
occasions  in  practice  when  it  is  of  interest  to  know  that  an  approximation 
to  the  solution  of  the  Dirichlet  Problem  can  be  trusted  on  the  boundary. 

In  the  discussion  immediately  above  it  was  brought  out  that  even  with 
the  smoothest  conceivable  C  (a  circle)  some  smoothness  conditions  are 
still  needed  on  u  beyond  continuity  to  guarantee  con’/ergence  on  the  boundary. 
It  is  of  some  theoretical  interest  to  note  that  If  we  are  willing  to 
consider  the  smoothest  conceivable  u,  then  C  can  be  an  entirely  arbitrary 
Jordan  curve  and  convergence  will  take  place  on  C.  We  now  express  this 
result  formally.  The  mapping  function  which  appears  in  the  statement 

is  the  one  used  above  to  define  curves  of  Type  W, 

Theorem  3.5  .  l£t  D  be  a  region  V'ouncied  by  an  arbitrary  Jordan  curve  C. 
let  the  function  U  be  harmonic  on  the  closed  set  D  ♦  C ,  Let  denote 
the  image  in  the  z-planc  of  the  circle  [w  |  -  R  >  1  under  the  conformal 

map  given  by  z  ■  ij  (w) ,  Then  (a)  there  exists  a  largest  value  of  R, 

say  OD  ,  such  that  U  with  its  possible  harmonic  extensions  is  single- 
valued  and  harmonic  at  every  point  interior  to  j  (b)  there  exists  a 
sequence  of  point  sets  S^,  S2,  ...  on  C  such  that  for  any  R,  1  <  R  <  'P  , 

where  is  the  harmonic  polynomial  of  degree  at  most  n  foiuid  by  inteipolar- 
tion  to  U  in  the  podnts  S  . 


In  other  terms,  the  degree  of  convergence  of  the  sequence  of  harmonic 
interpolation  polynomials  on  D  +  C  is  Oil/f^),  for  juiy  R,  1  ■  R  <  "P , 
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Statement  (a)  in  the  conclusion  of  the  theorem  is  a  consequence  of 
the  monotonic  character  of  the  level  curves  C^,  For  details  the  reader 
is  referred  to  [17,  chap*  IV], 

For  statement  (b),  we  make  use  of  a  theorem  of  Walsh  [llt][l6]  which 
states  that  there  exists  a  sequence  of  harmonic  polynomials  h^,  hg,  *••, 
hjj,  •••,  of  respective  degrees  at  most  1,  2,  n,  such  that 

(3.9)  |U(2)  -  h  (z)|  <  4 

n  - 

for  z  C  and  for  any  R,  0<R<-f,  where  M  is  independent  of  z . 

We  let  be  the  point  set  referred  to  in  Theorem  3.2.  Substituting 
into  (3.2)  we  obtain 

lH^(Uj  z)  -  U(z)  I  <  h  —  ,  zeD+C. 

R 


Taking  the  n-th  root  of  bothi  sicies  and  tlien  passing  to  the  limit,  we 
obtain  the  conclusion  of  Theorem  3.?,  since 

Walsh's  harmonic  approximation  result  embodied  in  (3.9)  is  valid 
when  D+C  is  replaced  by  any  closed  bounded  point  set  P  whose  complement 
is  connected  and  regular  in  the  sense  that  it  possesses  a  Green 's  function 
with  pole  at  infinity'.  For  details  see  [16]  and  [17,  chap,  b].  Theorem 
3.5  can  be  generalized  correspondingly,  provided  that  a  hypothesis  is 
inserted  to  the  effect  that  for  each  n,  P  has  a  subset  satisfying 
the  n-s  condition. 


By  using  interp-oiation  techniques,  Shen  [17,  pp.  173-17ll]  constructed 
a  proof  of  the  existence  of  sequences  of  complex  polynomials  converging 
maximally  (the  terminoiogy  is  that  of  Walsh)  to  a  given  complex  function 
f  analytic  on  a  general  closed  bounded  point  set  whose  complement  is 
connected  and  regular,  Shen's  polynomiais  were  defined  by  interpolation 
to  f  in  points  with  extremal  properties  similar  to  the  extremal  properties 
considered  in  Theorem  3.2  above.  It  would  be  of  interest  to  see  if  with 
the  aid  of  Theorem  3.2  an  analogous  existence  proof  could  be  given  for 
harmonic  polynomials  and  harmonic  functions.  Theorem  3*5  and  the  other 
theorems  of  this  section  of  course  are  not  existence  theorems  in  the 
basic  sense  since  they  presuppose  the  existence  of  some  sort  of  harmonic 
polynomial  approximations. 
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U •  Harmonic  polynomial  interpolation  In  transforms  of  the  roota  of 
unity.  We  now  consider  the  following  infinite  sequence  of  point  sets 
lying  on  a  Jordan  curve  C  in  the  *-plane8 

Sj  -  {  ll)  (wj),  ij)  (w^),  ijl  (1)  ), 

Sg  -  Mw^),  lll(w^),  l|j(w^),  MD), 

0 

.  *<’‘2^1)- 

• 

here  ''2n+l  "  e^’'^/^2n+l)  z  ■  ij)  (w)  is  the  mapping  function  which 

in  Introduced  in  Section  3  to  define  curves  of  Type  W,  (That  is,  t|J  is 
r.alytic  and  univalent  on  |w|  >  1,  continuous  on  |w|  >  1,  and  maps 
n'  1  >  1  onto  the  region  exterior  to  C  so  that  the  points  at  infinity 
irrespond.)  It  is  known  that  this  sequence  of  interpolation  points  is 
^  fundamental  importance  in  complex  polynomial  interpolation  to  boundary 
•hif's  (see  [U]),  so  it  is  natural  to  study  it  in  connection  with 
rmonic  polynomial  interpolation. 

Let  u  be  a  continuous  function  on  C,  We  choose  an  infinite  sequence 
'o?!se  polynomials  p^,  P2,  p^^,  of  respective  degrees  1,  2, 

;  we  assume  that  satisfies  the  n-s  condition  of  Section  2, 

I  we  set  up  the  harmonic  interpolation  polynomial  in  the  form  (2.U)  or 

Ul)  : 

2n+l  . 

Ii.l)  H^(u,  z)  .  Ut  l|l 

.’here 

■lu?)  B^^(z)  -  .  ^£jb<5>  pj(z)  ppr] 

■  •’ko’  *  f 

!  .notion  is  the  unique  harmonic  polynomial  of  degree  at  most  n 

I  )c 

wdshes  at  all  points  of  S^  except  ij)  equals  unity. 

'  hic  above  representations  of  and  arc  not  quite  as  general  as 
■  f ,h t  be,  in  that  it  wou3.d  be  theoretically  possible  to  choose  a  new 


set  of  base  polynomials  each  time  n  is  changed.  We  might  then  use  the 
notation  p^^,  p^g*  •••»  Pnn  generality  is  not 

needed  in  what  follows. 

We  shall  now  calculate  Hj^(u|  *)  explicitly  for  the  special  cases 
in  which  C  is  a  circle  and  C  is  an  ellipse. 

For  the  circle  |e|  ■  R,  the  mapping  function  becomes  b  -  l|)  (w)  ■  Rw 

2  2  n  a_n 

We  choose  as  the  base  polynomials  Pj^  ■  b/R,  P2  ■  b  /R  ,  ...)  p^^  ■  /tr ^ 
The  coefficients  of  ^re  the  solution  of  the  following  system  of 

linear  equations  in  the  unknowns  ^kj^  " 


(U.3) 


>1  J 


^  w 

"2n+l 


h  -  1,  2,  ...»  2n  1, 

whor«  *0,  h  y  k,  6j^j^  •  1.  The  existence  of  a  unique  solution  to 

(I1.3)  will  be  a  by-product  of  our  method  of  finding  the  solution.  We 

shall  need  the  following  easily  proved  facts  about  the  (2n+l)— th  roots 
h 

of  unity,  h  ■  1,  2,  ...,  2n+li 

(a)  '^2n+l  ^2n+l  “  integers  J  and  h. 

0,  3^0(mod2n+l) 

2n+  1,  J  =0(mod  2n+  l). 


2n+l 

I  w  "J 
.  2n+] 

h“l 


We  now  transform  the  system  (U»3)  into  another  much  simpler  system 
as  follows.  First  we  add  ail  the  equations  together.  Duly  referring  to  . 
(b)  above,  we  find  that  the  sum  equation  is  merely  the  equation 
(2n+  l)b^  -  1,  Then  for  each  J,  j  ■  1,  2,  n  we  multiply  the  first 

equation  by  the  second  by  (w2n+i^  »  •••>  ^^e  (2n  +  l)-thby 

and  add  the  resulting  equations.  We  find  by  referring  to 

(a)  and  (b)  that  the  sum  equation  reduces  to  (2n  ♦  l)bj  ■  ^2n+l* 

s.™  procedure  with  Inetead  of  (2n.  1)B,  -  .  J",. 
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In  matrix  language,  this  process  simply  amounts  to  replacing  the 

linear  system  (U*3),  which  has  a  matrix  we  shall  here  denote  by  W,  by 

a  new  system  of  linear  equations  In  which  the  matrix  Is  W  W  and  the 

right  side  Is  the  k— th  column  of  W  ,  where  W  Is  the  conjugate  transpose 

of  W,  The  matrix  produce  W  W  Is  a  scalar  matrix  with  diagonal  elements 

— ’1/2 

all  equal  to  2n+l.  The  matrix  (2n+l)  '  W  is  unitary.  The  absolute 

value  of  the  determinant  of  W  may  be  calculated  as  follows! 

det  W*  W  -  HiFW  det  W  -  ldetW|^  -  (2n  + 

(U.U)  |det  W|  -  (2n  + -  |detW*|. 

This  result  is  of  some  interest  in  connection  with  the  proofs  of 
Theorems  3.2  and  3.3 »  In  which  the  determinant  of  the  matrix  for  the 

"  I  ^ 

analogous  general  case  is  maximized. 


The  particular  consequence  of  the  above  discussion  needed  here  is 
that  W  is  non-singular,  so  (U.3)  has  a  unique  solution  and  it  is  moreover 
the  same  as  the  solution  to  the  new  system  having  the  matrix  W  W. 


The  solution  of  {U.3)  is  now  obvious, 
obtain 


'nlc'*' 


Substituting  into  (U.2)  we 


f  —  k  \  4 

n 

1  +  S 

(^) 

_ 

2n  +  1 


-By  summing  the’ geometric  series,  this  can  be  written  as 


*  Suppose  that  in  (U.3)  we  replace  h  -  1,  ...,  2n+l  by  any  set 

of  2n  + 1  points  on  the  unit  ■  circle.  The  Hadamard  determinant 

inequality  yields  the  middle  member  of  (U.U)  as  an  upper  bound  for 

|det  W I  no  matter  how  these  points  are  chosen.  Therefore  the  choice  of 
points  as  the  roots  of  unity  has  maximized  (det  W I . 
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If  z  lies  on  some  closed  circular  disk  (z|  <  R*  <  R,  the  first  term 
in  braces  is  positive  but  less  than  unity  and  the  second  term  approaches  aero 
with  n  at  a  geometric  rate  of  convergence.  Therefore  ■  0(l/(2n+l)) 

■  0(l/n)  uniformly  ink,  k«l,  ...,n,  and  uniformly  for  (z(  <  R*. 

It  follows  from  the  Corollary  of  Theorem  3«1  that  with  (U,5)  substituted 
into  (U.l),  jj^iy  ^  exists  uniformly  on  |z|  <  R*  <  R,  and 

provides  there  a  solution  of  the  Dirich  let  problem  for  u  and  the  region 

|a  I  <  R. 

It  is  easy  to  hyjKiss  Theorem  3«1  io  the  present  case  in  demonstrating 
the  convergence.  We  write  *)  i’^  i*''®  form 


(b.6) 


1 

7^ 


2n4-l 

L 

k-l 


u(R  w 


^  ) 
2n+l^ 


r2- 


2ri+l 


2n 

2n  +  1 


2R 


The  first  summation  term  is  for  [a |  /  R  a  Riemann  sum  approximating 
the  Poisson  integral  fl,  p.  180] 

I.i2 

U(z) 


1  ,  iSx 

^  /  u(R  e  )  - 

0  R 


-  z 


which  represents  the  solution  of  the  Dirichlet  Problem  for  u  and  the 
disk  |zl  <  R.  It  is  an  easy  exercise  in  analysis  to  show  that  the 
second  summation  term  approaches  zero  for  |z|  <  R,  so  again  we  have 
lim  H  (u;  z)  =  U(7,),  |z|  <  R. 

This  argument  establishes  convergence  on  |zl  <  R  when  u  is  merely 
Riemann  integrable  and  not  neressarlly  continuous . 

On  the  circle  [z|  ”  R  itself,  a  calculation  shows  that 


(U.7)  nj^^(R  e^®) 


sin(n  ^  •i)(e  - 


(Jn  .  ])  .In 


and  when  thin  is  substituted  into  H^(u;  z) ,  wo  obt,ain  the  trigonometric 
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J  A 

sum  of  degree  at  most  n  found  by  interpolation  to  u(R  e  )  in  the  points 

S.  -  2irk/(2n+l),  k  ■  1,  2n+  1.  The  identification  of  H  (uj  z) 

ic  n 

on  the  circle  with  a  trigonometric  sum  was  previously  mentioned  in 
Section  3»  As  stated  in  Section  3»  it  is  known  [6,  p,  120]  that 

T  (R  e^®)  -  ElB.(Re^®)l  -  O(logn). 

n  nic 

From  (3 *2)  and  the  facts  given  in  the  discussion  following  Theorem  3*1 
it  follows  that  if  u  satisfies  a  Upschltz  condition  of  order  o<,  0  <  o<  <  1 
on  |zl  -  R,  then  H^(z)  —  U(z)  -  0(log  n/n*^)  uniformly  on  |z|  <  R, 
so  convergence  takes  place  uniformly  on  the  closed  disk.* 


We  turn  to  the  case  in  which  C  is  an  ellipse,  which  we  shall  call  E. 

2  2  2  2 

If  we  let  z  •  X  ♦  ly,  the  familiar  equation  (x  /a  )  +  (y  A  )  ■ 

a  >  b  >0,  may  come  first  to  mind,  but  here  it  is  more  convenient  to 

define  positive  constants  R  and  c  respectively  by  R  •  ^(a  +  t))/^a  —  1}), 

/  2  ^ 

c  ■  /a  -  b  and  write  the  equation  of  E  in  the  form 


(lu8) 


c^(R 


c‘(R-^) 


1,  c  >  0,  R  >  1, 


The  mapping  function  for  E  is  given  by  [l,  pp,  76-77] 


(U.9) 


c(R.  *  i). 


w|  >  1, 


This  mapping  function  is  analytic  and  univalent  for  (w{  >  l/R  and 
gives  a  conformal  map  of  this  region  onto  the  region  in  the  z-plane 
which  consists  of  the  entire  plane  minus  a  cut  on  the  real  axis  from  -c 
to  +c.  In  this  map,  the  exterior  of  any  circle  (w|  ■  "f  >  lA  is  mapped 
onto  the  exterior  of  a  certain  ellipse  in  the  z-plane  having  an 
equation  like  (h.8),  but  with  R  replaced  by  'fR.  It  is  then  easily  seen 


*  The  explicit  formula  (lj.5)  may  have  been  first  published  only  in  I960  [U], 
but  both  it  and  the  convergence  results  noted  above  were  certainly  known 
to  Walsh  in  the  late  1920 *s.  He  gave  the  convergence  theorem  to  the 
author  as  a  student  exercise  in  1932,  and  in  [13]  and  [lij]  he  pointed  out 
the  relation  between  harmonic  polynomial  interpolation  on  the  circle  and 
trigonometric  interpolation. 
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that  if  'f  <  1,  then  must  lie  inside  E,  and  that  given  any  point 
inside  E,  for  some  sufficiently  near  to  1  this  point  also  lies  inside  E^, 

Henceforth  in  the  discussion  we  restrict  w  to  the  region  Iwj  >  1/K, 

In  this  region  the  transformation  (U.9)  carries  Rw  +  l/(Rw)  into  the 

polynomial  z/c  and  (Rw)  +  l/(ftw)  into  the  polynomial  (z/c)  -  2,  We 

shall  abbreviate  a  statement  of  this  type  by  saying,  for  example,  that 

(Rw)^  +  l/(Rw)^  "is"  a  polynomial  in  z.  Now  if  it  is  true  that 

(Rw)®  +  l/(Rw)®  is  a  polynomial  in  z  of  degree  m  for  m  ■  N  —  1  and  ■  ■  N, 

then 

(Rw)^  ^  ■  [(Rw)^  — i-n)(Rw  ♦  ik)  -  [(Rw)**”^  +  — -Tijr] 

(Rw)”*^  (Rw)“  ^  (Rw)^^ 

is  a  polynomial  in  z  of  degree  N+  1,  The  hypothesis  is  true  for  N  ■  2, 
so  we  have  established  that  p^  -  (Rw)®  +  l/(Rw)®  is  a  polynomial  of  degree 
m  in  z  for  m»  1,  2,  ...  . 

We  take  the  polynomials  p^  as  our  base  polynomials  in  constructing 
(U.l)  and  (U.2)  for  the  ellipse.  The  points  in  are  of  course  now  the 
points  +  l/(Rw2n>i)J»  2,  ...,  2n  ♦  1, 

The  coefficients  bj,  of  B^)j(®)  are  the  solution  of  the  system  of 
linear  equations 


h  •  2^  2n + 

To  solve  the  system,  we  premultiply  the  matrix  of  the  system  and  also  the 
vector  on  the  right  hand  side  by  the  matrix  W*  introduced  above  in  the 
circle  case.  This  time  we  arrive  at  the  following  system  of  linear  equations: 

(2n  l)b  -  1 

0 

(2n  +  1)[R^  ^  J  2,  ...,  n 

(2n  +  1)[-^  bj  J  “  1»  2,  ...,  n. 


(U.IO) 


n 

♦  Zb 


>1 


!<■>  • 


(R  W  )^ 
^  2n+l^ 


n 

♦  Zb. 
J-1  ^ 


1 

~=ir 


«2‘n^l)' 
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k 

’2n+l* 


where  Direct  calculation  shows  that  the  determinant  of  the 


transformed  system  has  the  value  (2n  +  it  (R^*^  -  RT^^).  If 

J*  ^ 

A  denotes  the  matrix  of  the  original  system  (U.IO),  then  det  W  det  A 
has  this  value,  so  by  (U.b), 

|det  A I  -  (2n  ♦  i)(2n+l)/2  ^  ^j^2j  _  1  ^  ^  ^  ^ 

j-1  R 


The  value  of  the  determinant  is  of  some  interest  in  connection  with  finding 
interpolation  points  with  extremal  properties  in  connection  with  Theorems 
3.2  and  3.3. 

The  discussion  of  the  determinant  of  course  implies  that  (U.IO) 
has  a  unique  solution  which  can  be  found  from  the  transformed  system. 
Accordingly  we  find  that  the  solution  is  in  part 

1 

2n  ♦  1 


(rO)J - L_ 

(RU>)J 

^  J  •  1,  n, 

(2n+l)(R^J  -  ) 


and  the  complex  conjugates  of  the  numbers  bj  complete  the  solution. 
Substituting  into  (U.2),  we  obtain 

(RW  - 


(Uai)  (2n»l)B  .(i)  -  1  ♦  5?2  £  — - t(lh.)J  t -i-, 

>1  -  4,  («v^ 


where  as  above  and  |w|  >  1^, 


To  study  the  order  of  magnitude  of  6.(2)  as  n  becomes  infinite, 

4  0  r%  TiK 

we  let  w  ■  '^e^  ,  y  >  1/R,  and  rewrite  (U.ll)  in  the  form 


(li.l2)  (2n  +  l)B^Jz)  -  1  +  R.2  I 

j-1 


_  _! 
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From  this  it  follows  that 


1  +  2  Z 

J-1 


(U.13)  |B„^(*)|< 


-  -i. 


2n  +  1 


Now  let  a  lie  oh  the  locus  which  is  the  map  of  |w|  ■  "f^, 
where  (l/R)  <  'P  <  1.  It  is  possible  to  remove  the  absolute  value 
Indications  in  (U.13),  and  when  this  is  done  the  right  side  becomes 


n 

1  +  2  Z 

>I 


2n  +  1 


(B 'F  )  J  +  .1 

n 

1  +  2  Z 

J-1 

inf)i 

R'^  +  i 

RJ 

2n  +  1 


n 

1  +  2  Z 


(Rf 


2n  +  1 


1^2  Z  (f J 

j-1  R*’ 


2n  +  1 


The  two  geometric  progressions  in  the  n\imerator  of  the  last  member 
converge,  so  we  have  shown  that  “  0(l/(2n+l))  -  0(l/n)  uniformly 

for  z  on  .  But  since  is  a  harmonic  polynomial,  it  follows  by 
the  maximum  principle  that  -  0(l/n)  uniformly  for  z  on  and 
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Interior  to  Ef  •  It  will  be  recalled  that  nan  be  adjusted  so  as  to 
contain  any  closed  point  set  Interior  to  E»  Therefore  according  to  the 
Corollary  of  Theorem  3.1*  when  H^(uj  z)  In  (U.l)  Is  constructed  with 
®nk^*^  as  given  by  (U.ll)  or  (U.12),  ^(u|  z)  exists  uniformly 

on  any  closed  point  set  of  the  Interior  of  the  ellipse  (U»8)  and  provides 
there  the  solution  of  the  Dlrlchlot  problem  for  the  boundary  data  u. 

Going  back  to  (U.12),  It  can  be  seen  that  If  •  1,  so  that  z  lies 
on  the  ellipse  E,  then  reduces  to 


1.21  cos  1(8  - 


2n  +  1 


sin  (n 

(2n  +  l)  sin  ^(9  - 

This  Is  the  same  as  the  formula  (U.7)  which  we  obtained  In  the  circle  case, 
and  once  again  a  Llpschitz  condition  on  u  will  Insure  convergence  In  the 
closed  region  under  consideration. 

The  identification  of  z)  for  z  on  the  ellipse  with  trigonometric 

interpolation  was  acconiplished  by  Walsh  [l8]  who  then  pointed  out  the 
convergence  implication  under  suitable  smoothness  conditions  on  u.  The 
present  discussion  eliminates  any  requirement  on  u  beyond  mere  continuity 
for  convergence  interior  to  the  ellipse. 

For  convenience  in  reference  we  summarize  the  above  results  formally. 
Theorem  U,1  .  Let  C  be  either  the  circle  |z|  ■  R  or  the  ellipse 

-  1,  c  >  0,  R  >  1, 

Let  the  function  u  be  continuous  on  C,  Let  *)  be  the  harmonic 

polynomial  of  degree  at  most  n  which  coincides  with  u  in  the  points 
®nk  "  ^’*2n+l^»  k  =  1,  2,  2n+  1,  where  Is  a  (2n  +  l)-^ 

root  of  unity  and  l|l  gives  the  conformal  map  of  |w|  >  1  onto  the 
exterior  of  C  so  that  the  points  at  infinity  correspond.  The  polynomial 
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exists  and  Is  vuiiquely  determined,.  For  the  circle  H^(u;  z)  Is  given 
by  (U«6)  for  z  not  on  C,  and  for  the  ellipse. 


2n+l 

Z 

k-1 


n 


(R  OJ )^  - 


1*%2  Z 
J-1 


jr 


- [(Rw)J  +  -JU] 

-  -L  (Rk)J 


2n 


$ 


c*j  ■  w, 


2n-*-l* 


c(Rw 


1 

1 


<  w 


In  either  case,  lim  h  (uj  %)  exists  uniformly  on  any  closed  subset  of 
-  n~»  00  n  *  - “ - 

the  Interior  of  C  and  provides  there  the  solution  of  the  Dirichlet  problem 
for  the  data  u.  Also  In  either  case,  for  z  on  C  and  with  z  ■  ijl  (e^^), 

H  (u;  z)  reduces  to  the  tri^'onometric  sum  of  order  at  most  n  which 

n  - - - jg - -  - 

Interpolates  to  u(l)  (e  ))  in  the  points  Oj^  "  2nk/(2n-*-  1),  k  ■  1, 

2n +  1.  If  u  satisfies  a  Lipschitz  condition  of  order  o<,  then  JiSooV”'  •> 
converges  uniformly  on  and  inside  C, 

It  might  be  of  interest  to  devise  a  proof  of  convergence  in  the 
ellipse  case  from  (h.lli)  directly  in  such  a  way  that  the  continuity  condition 
on  u  could  be  relaxed  to  an  integrability  condition  of  some  sort. 


We  conclude  with  a  technical  note  looking  toward  the  generalization 
of  the  above  theorem  to  more  or  less  arbitrary  Jordan  curves  C.  In 
constructing  for  the  ellipse,  we  could  just  as  well  have  used 

p*  -  c^((Rw)J  ^  l/(Rw)J),  J  >=  1,  ..  . ,  n,  as  our  base  polynomla]r  in  z« 

J  1 

After  cancelling  a  c  which  would  have  appeared  in  the  numerator  and 
denominator  of  each  term  in  the  summation  in  (U.ll),  the  formulas  (U.ll), 
(li*12),  and  (U.13)  would  have  looked  just  the  same  as  they  do  now.  These 
polynomials  pj  are  the  so-called  Faber  polynomials  for  the  region 
exterior  to  the  ellipse  E.  In  the  case  of  the  circle  |z|  ■  R,  the  Faber 


*  See  [15,  pp.  32-33]  for  a  quick  summary  and  a  number  of  primary  references. 

There  appears  to  be  a  slight  misprint  in  [1^]  in  the  recursion  formula 

at  the  bottom  of  page  32,  where  the  term  a  should  be  (n  l)a  . 

°  ’  n  n 
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polynomials  are  z,  z^,  z”,  which  are  equivalent  for  present 

purposes  to  the  base  polynomials  which  we  used.  Our  method  for  finding 

H^(u}  z)  in  each  case  thus  consisted  in  expressing  B^j^(z)  as  the  real 

part  of  a  linear  combination  of  the  first  n  Faber  polynomials,  and  then 

adjusting  the  coefficients  of  the  linear  combination  so  that  the  real 

part  would  vanish  at  2n  of  the  interpolation  points  and  equal  unity  at 

the  remaining  point.  The  orthogonality  properties  of  the  values  of  the 
2  —  -2 

monomials  1,  w,  w  ,  w,  w  ,  in  the  roots  of  unity  were  used 

in  the  adjustment  process. 

In  the  general  case,  the  mapping  function  z  ■  (w)  can  be  chosen 

so  as  to  have  a  Laurent  development 

(U,15  z  -  rw  +  a^  +  —  +  +  ...  ,  r  >  0,  |w|  >  1, 

w 


where  r  is  the  transfinite  diameter  or  capacity  of  the  curve  C,  The 
Faber  polynomials  for  C  are  defined  by  stating  that  Pjj(z)  must 

be  of  degree  n,  must  have  unity  for  the  coefficient  of  z*',  and  must  be 
such  that  the  Laurent  development  of  as  a  function  of  w  is  of  the  form 


Pn(.) 


/  \  n 
(rw) 


^n2 

“T 


+ 


•  •  • 


|w|  >  1. 


In  the  case  of  the  ellipse,  all  of  the  coefficients  g^j  vanish  except 
g^,  and  thia  paved  the  way  to  obtaining  a  simple  explicit  formula  for 
the  required  linear  combination  of  1,  Pj^(z),  PgCz),  ...,  P^(z),  P^(z), 
T^{z)  f  ...,  l*^(z)  .  In  more  general  cases  one  cannot  hop)e  to  obtain 
manageable  explicit  formulas  in  terms  of  the  coefficients  in  (U.l5),  but 
it  may  be  possible  to  establish  the  asymptotic  properties  of 
by  using  known  asymptotic  properties  of  the  polynomials  [15,  p,  33] 
together  with  the  orthogonality  of  the  sequence  on  the  unit  circle. 
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